Wigner function for squeezed negative binomial state and evolution of density operator for amplitude decay
Lv Heng-Yun1, Wang Ji-Suo2, Zhang Xiao-Yan1, Wu Meng-Yan1, Liang Bao-Long1, Meng Xiang-Guo1, †
Shandong Provincial Key Laboratory of Optical Communication Science and Technology, School of Physical Science and Information Engineering, Liaocheng University, Liaocheng 252059, China
Shandong Provincial Key Laboratory of Laser Polarization and Information Technology, College of Physics and Engineering, Qufu Normal University, Qufu 273165, China

 

† Corresponding author. E-mail: mengxiangguo1978@sina.com

Abstract

Using the thermal-entangled state representation and the operator-ordering method, we investigate Wigner function (WF) for the squeezed negative binomial state (SNBS) and the analytical evolution law of density operator in the amplitude decay channel. The results show that the analytical WF is related to the square of the module of single-variable Hermite polynomials, which leads to a new two-variable special function and its generating function, and the parameters s and γ play opposite roles in the WF distributions. Besides, after undergoing this channel, the initial pure SNBS evolves into a new mixed state related to two operator Hermite polynomials within normal ordering, and fully loses its nonclassicality and decays to vacuum at long decay time.

1. Introduction

Squeezed states and their practical applications have been widely explored during the past several decades (see, e.g., the review articles[1,2] and references therein). This is mainly because the noises of squeezed states may be below the vacuum or ground-state noise level in a certain observable, which can be useful for interferometric detection and gravitational waves by means of improving the interferometric and spectroscopic measurements. Especially, two-mode squeezed states served as a new entangled resource is crucial in realizing some vital quantum information processing, such as quantum teleportation,[3,4] quantum computation,[5] and cryptography,[6] and can be readily generated by using experimentally controllable and available sets including phase shifters, beam splitters, and squeezers.[7] For another side, any system in the nature is not absolutely isolated, quantum noise usually occurs and causes the decoherence of system when this system is in the interaction with its surrounding environment. Hence, the environment-induced decoherence aspects in open squeezed optical fields have received great attention from researchers. For instance, the effect of amplitude decay, phase damping, and thermal noise on the decoherence of squeezed states has been analytically investigated in the pioneering works.[810]

In probability theory and statistics, the negative binomial distribution reads , where and s is a non-negative integer. Corresponding to this definition, there exists a pure state with a negative binomial form of its photon-number probability distribution, which is called the negative binomial state (NBS),[11] that is

where is the Fock state obtained by repeatedly acting the photon-creation operator on the vacuum state in the Fock space. Using the identity , we can rewrite as
Indeed, the NBS is an intermediate state between a coherent state ( and ) and a quasi-thermal state (s = 0 ). The reason for designating the state for s = 0 as a quasi-thermal state is that its photon number distribution is the same as that of a mixed thermal state, which is easily proved to be true using the compact form , regarded as a multi-photon-subtracted thermal state via subtracting s photons from a single-mode thermal field with the mean number of photons .[12] This means, through subtracting s photons, the photon number of shows a negative binomial distribution. For the NBS and its derived states, various representations, properties, and quasiprobability distributions have been studied in some detail and their generation schemes in a number of nonlinear processes have been demonstrated.[1315] Especially lately, the applications of the NBS in quantum optics and its time-evolution in some typical decoherence channels including the amplitude decay channel (ADC),[16] the diffusion process,[17] and the laser channel[18] have been investigated.

As a practical derived state, a squeezed negative binomial state (SNBS) was first introduced by Joshi and Lawande, and its second-order correlation function, Wigner function, Q functions, and photon number distribution have been investigated, and one result that needs to be noted is that the SNBS can be helpful for the transfer of quantum phase-space distributions from squeezed coherent states to squeezed number states.[19] In this paper, we shall introduce a new formula of SNBS, which is defined as an NBS having undergone a two-mode squeezing transform, hence its density operator reads

where is the single-mode squeezing operator with the squeezing r.[2022] To go a step further, we focus on the theoretical investigation of Wigner function (WF) for SNBS and the evolution behavior of density operator in ADC by fully using the operator-ordering method[2326] and the thermal entangled state representation (TESR).[2730] Our motivations for focusing on this topic are explicated as follows. First, the operator-ordering method can help us prove that the WF for SNBS is expressed as the finite-sum representation determined by the non-negative integer s, rather than the infinite-sum form as shown in Ref. [19], which can more accurately describe the nonclassicality of SNBS in phase space (see Section 3 for details). Second, according to our best knowledge, the investigation of the evolution of SNBS in ADC by deriving the analytical results has not been previously discussed because of the complicated calculations resulting from the squeezing operation in SNBS, however TESR can help one to greatly simplify it and derive an accurate and concise result related to two operator Hermite polynomial functions.

2. Operator-ordering products of

To analyze the nonclassicality of SNBS from the aspects of phase space and decoherence evolution, we first use the Bogolubov transformation caused by the squeezed operator to change SNBS into

where , a weighted quantum superposition of squeezed number states . Next, we derive the normal and antinormal ordering products of . For this purpose, using the completeness relation of the coherent state to derive the normal ordering of the operator , that is
where the marked symbol refers to normal ordering.[31] Further, using the integral formula (17) in Ref. [32] to perform the integration (5) over α within normal ordering, and operating this integral result on the vacuum state and using the generating function of n-order Hermite polynomial , we have the compact form of the squeezed number state , that is
which indicates that a squeezed number state can be viewed as a Hermite polynomial excitation on the squeezed vacuum state. Noting the normal ordering product of the vacuum projector, i.e., , we directly obtain the normal ordering of the density operator , i.e.,
Using the summation formula concerning two single-variable Hermite polynomials in Eq. (10) in Ref. [33], thus the normal ordering expansion of is obtained as the compact exponential form
where we have set , , and . Using a method similar to that used to obtain the result (20) in Ref. [33], we can also prove that the normal ordering operator is a bivariate distribution operator in phase space. Further, noting the operator identity , equation (8) can be rewritten as
which shows that the operator is indeed a squeezed thermal state by comparing with the result (11) in Ref. [33].

Further, substituting the normal ordering of in Eq. (8) into the antinormal ordering formula[34] and using the mathematical integral formula,[35] we can obtain the antinormal ordering of as the following form

where denotes antinormal ordering[36] and we have set , , and .

3. Wigner function for SNBS

The Wigner distribution function is a useful tool which is used to investigate the nonclassical optical fields from the view point of phase space. In this section, we shall investigate the nonclassicality of the SNBS by analytically and numerically studying its phase-space WF distributions. For the given density operator , its WF definition is given by , where is the coherent state representation of single-mode Wigner operator, i.e., .[37] Considering that are coherent states,[38,39] hence we first introduce the normal operator identity of the operator ,[10] that is

and use the antinormal ordering of in Eq. (10), thus the state reads
where and . Thus, inserting Eq. (12) into the WF definition and using the completeness relation of the coherent state and the generating function of , i.e.,
we therefore have

Further, using the integral formula[40,41] to complete the two integrations over β and z in the whole phase space, the WF can be calculated as

where we have set
with Further, using the exponential expansion of cross product term , it then follows that
Noting the generation function of in Eq. (13) and its l-order differential relation of , we finally obtain
which is the analytical WF of the SNBS related to the square of the module of -th single-variable Hermite polynomial. Obviously, when s = 0, becomes the WF for squeezed quasi-thermal field, i.e.,

Mathematically, it is well known that two-variable Hermite polynomials can be introduced either by its generating function

or by its power-series expansion
where represents the transfer of the amplitude of number state in the time evolution process of a forced Harmonic oscillator,[42] and can be useful for studying light propagation modes in quadratic graded-index media.[43] Hence, a comparison of the right sides of Eqs. (15) and (18) directly leads to the result that

If we replace , , and in Eq. (22), we therefore have

where
is a new two-variable special function whose generating function is , a new and interesting result. Thus, the differential relation of the special function can be obtained as
which is the same form as the well-known differential relation of .[29] Further, using the special function , we find that the right side of Eq. (22) can be represented as , thus the analytical WF for SNBS can be simplified as

To clearly observe the nonclassicality of SNBS, we plot its WF in (Re ) phase space for a fixed squeezing r and different values of the parameters s and γ as shown in Fig. 1. Obviously, the SNBS always shows squeezing effect in one of the quadratures for as a feature of nonclassicality of this state. For (e.g., s = 100) and (or s = 0), the WF becomes the standard Gaussian distribution of squeezed coherent state (or squeezed quasi-thermal state) [see Figs. 1(a) and 1(d)], which is completely accordant with the analytical result outlined in Ref. [19]. Comparing Figs. 1(a), 1(b), and 1(c), we find that, with the increase of s, the WF gradually shows an inverted peak in the center of phase space, however it slowly loses the inverted peak as γ increases and finally becomes Gaussian when γ tends to one. Moreover, the larger value of s can lead to the smaller peak value of WF, but the variation of the peak value with γ is just a fully contrary trend. This result can easily be understood via analyzing the result in Eq. (2) since more and smaller number states are annihilated and the probabilities of finding number states in this field becomes small with increasing s (especially at center part in the phase space), however the larger γ causes the higher weights of the existing number states. Besides, for the limited case ( ), the WF distribution seems basically unchanged with the increase of s, which is because the value of parameter s has a negligible effect on the weights of the existing number states.

Fig. 1. The WF distributions for SNBS with a fixed squeezing r=0.4 and the different values of the parameters : (a) (0,0.5), (b) (5,0.5), (c) (5,0.9) and, (d) (100,1).
4. Evolution of SNBS in ADC

For any open system, noise always arises in fundamental dynamics processes. For instance, the purely dissipative noise in ADC leads to the nonclassicality deterioration of the system. In this section, we shall analytically obtain the evolution formula of SNBS in ADC and investigate how does the amplitude dissipation affect the nonclassicality of SNBS.

4.1. Kraus operator solutions to master equation describing ADC

In the interaction picture, the density-operator evolution of the system in ADC can be characterized by[4447]

where κ is the decay rate of the system. Indeed, this process completely differs from the phase-decay process with a purely dephasing noise. To obtain the solutions of the above operator equation, we review the two-mode TESR describing quantum entanglement between system and heat reservoir, which is defined in enlarged Fock space as
where is a displacement operator with the complex variable η, the operators and are respectively the real and fictitious creation operators, which yield the commutative relations, i.e., and . Under the basis of the entangled state , there exist the following extremely useful interchange relations between the real and fictitious modes, i.e., , , and , which play a vital role of solving the master equation (27) since it can help convert the equation for the density operator in Eq. (27) into the evolution equation for the state vector , that is
whose formal solution is
Fruther, through calculating the scalar product , we find that it is proportional to the inner product between an exponentially decaying entangled state and the initial state vector , i.e., , and its proportionality coefficient is just a Gaussian decaying factor , . This point clearly reveals the decoherence behavior of the channel, namely, the amplitude decay. To take off the density operator from , we need to use the completeness relation , thus after strict derivations (for the detailed formula derivations see Ref. [27]), its infinitive Kraus operator-sum representation of as the solutions of the master equation (27) can be calculated as
where is the Kraus operator of . Hence, for a known initial state , the evolved density operator in ADC can be freely obtained from Eq. (31). Also, using the operator identity , we can prove the normalization condition .

4.2. Density-operator evolution

Noting the marshalling sequence of the operators , , and in Eq. (31), we need to introduce the anti-normal operator identity of the operator ,[32] that is

and use the normal ordering of in Eq. (8). Hence, the state reads
Substituting Eqs. (8) and (33) into Eq. (31) leads to the form
where we have set and . Inserting the completeness relations of coherent states and into Eq. (34),[38,39] and using the operator identity and the generating function of in Eq. (13),[48] we obtain
Using the integration technique within an operator-ordering product[4952] to perform the two integrals over α and β, then becomes
where can be represented as
since the normal ordering symbol within another normal ordering symbol can be canceled, i.e., , and the parameters h, h1, h2, and h3 respectively are Comparing the right sides of Eqs. (15) and (36), and using the new special function in Eq. (24), we immediately obtain the analytical evolution of SNBS in ADC, that is
which means that the pure SNBS evolves into a new mixed state in ADC, related to two operator Hermite polynomial functions whose variables are the linear combinations of creation and annihilation operators within normal ordering. Especially, when r = 0, thus g = 1, , and , equation (38) becomes the evolution of NBS in ADC, i.e.,
with , which shows that, after undergoing ADC, the resultant state is still an NBS, but the parameter γ decays to and . It is in agreement with the result in Ref. [16]. For the case of s = 0, becomes
which is the evolution of squeezed quasi-thermal field in ADC. In addition, when the decay time , thus T = 0, h = 1, , , and , so equation (38) returns to SNBS in Eq. (3). For the limited decay time, i.e., , , thus SNBS finally decays to vacuum.

5. Conclusions

In summary, we have used TESR and the operator-ordering method to derive the analytical evolution law of SNBS in ADC. It is found that the analytical WF is the finite-dimensional sum over the square of the module of single-variable Hermite polynomials. Numerically, the WF for SNBS can exhibit a similar Gaussian distribution in both coherent and quasi-thermal limited cases, and the larger s can result in an inverted peak with the smaller peak value in the WF phase space but the larger γ plays an opposite role. For another case, after passing through ADC, the pure SNBS evolves into a new mixed state, whose normal ordering product of the density operator is related to two operator Hermite polynomials of the linear combinations of creation and annihilation operators, however the SNBS finally decays to vacuum as a result of amplitude decay when the decay time is long enough. From the above results, we have clearly found that, in the study of the environment-induced decoherence aspects, the TESR and the operator-ordering theory provide a privileged tool to investigate the analytical evolution of density operator and characterize destruction level of nonclassicality. Besides, through deriving the WF for SNBS, we have found a new two-variable special function and its generating function, with which the concise formulas of WF and density-operator evolution can be obtained.

Reference
1 Kimble H J 1992 Phys. Rep. 219 227 https://dx.doi.org/10.1016/0370-1573(92)90139-Q
2 Dodonov V V 2002 J. Opt. B Quantum Semiclass Opt. 4 R1 https://dx.doi.org/10.1088/1464-4266/4/1/201
3 Wang S Hou L L Chen X F Xu X F 2015 Phys. Rev.A 91 063832 https://dx.doi.org/10.1103/PhysRevA.91.063832
4 Sohma M Hirota O 2003 Phys. Rev.A 68 022303 https://dx.doi.org/10.1103/PhysRevA.68.022303
5 Barenco A Dutch D Ekert A Jozsa R 1995 Phys. Rev. Lett. 74 4083 https://dx.doi.org/10.1103/PhysRevLett.74.4083
6 Ekert A 1991 Phys. Rev. Lett. 67 661 https://dx.doi.org/10.1103/PhysRevLett.67.661
7 Mandel L Wolf E 1995 Optical Coherence and Quantum Optics Cambridge Cambridge University Press p. 649
8 Biswas A Agarwal G S 2007 Phys. Rev.A 75 032104 https://dx.doi.org/10.1103/PhysRevA.75.032104
9 Hu L Y Fan H Y 2008 J. Opt. Soc. Am.B 25 1955 https://dx.doi.org/10.1364/JOSAB.25.001955
10 Meng X G Wang Z Fan H Y Wang J S 2012 J. Opt. Soc. Am.B 29 3141 https://dx.doi.org/10.1364/JOSAB.29.003141
11 Berry M V 1984 Proc. R. Soc.A 392 45 https://dx.doi.org/10.1098/rspa.1984.0023
12 Fan H Y Zaidi H R Klauder J R 1987 Phys. Rev.D 35 1831 https://dx.doi.org/10.1103/PhysRevD.35.1831
13 Agarwal G S 1992 Phys. Rev.A 45 1787 https://dx.doi.org/10.1103/PhysRevA.45.1787
14 Matsuo K 1990 Phys. Rev.A 41 519 https://dx.doi.org/10.1103/PhysRevA.41.519
15 Wang X G Fu H C 2000 Int. J. Theor. Phys. 39 1437 https://dx.doi.org/10.1023/A:1003619724681
16 Chen F Fan H Y 2014 Chin. Phys.B 23 030304 https://dx.doi.org/10.1088/1674-1056/23/3/030304
17 Liu T K Wu P P Shan C J Liu J B Fan H Y 2015 Chin. Phys.B 24 090302 https://dx.doi.org/10.1088/1674-1056/24/9/090302
18 Da C Chen Q F Fan H Y 2014 Int. J. Theor. Phys. 53 4372 https://dx.doi.org/10.1007/s10773-014-2187-5
19 Joshi A Lawande S V 1991 J. Mod. Opt. 38 2009 https://dx.doi.org/10.1080/09500349114552111
20 Meng X G Wang Z Fan H Y Wang J S 2012 J. Opt. Soc. Am.B 29 1835 https://dx.doi.org/10.1364/JOSAB.29.001835
21 Xu Y J Meng X G 2013 Int. J. Theor. Phys. 52 4155 https://dx.doi.org/10.1007/s10773-013-1728-7
22 Wang Z Meng X G Fan H Y 2012 J. Opt. Soc. Am.B 29 397 https://dx.doi.org/10.1364/JOSAB.29.000397
23 Song J He R Yuan H Zhou J YFan H 2016 Chin. Phys. Lett. 33 110302 https://dx.doi.org/10.1088/0256-307X/33/11/110302
24 Fan H Y 2003 J. Opt. B: Quantum Semiclass. Opt. 5 R147 https://dx.doi.org/10.1088/1464-4266/5/4/201
25 Meng X G Fan H Y Wang J S 2018 Int. J. Theor. Phys. 57 1202 https://dx.doi.org/10.1007/s10773-017-3650-x
26 Zhou J Fan H Y Song J 2017 Chin. Phys.B 26 070301 https://dx.doi.org/10.1088/1674-1056/26/7/070301
27 Fan H Y Hu L Y 2008 Mod. Phys. Lett.B 22 2435 https://dx.doi.org/10.1142/S0217984908017072
28 Fan H Y 2008 Opt. Commun. 281 2023 https://dx.doi.org/10.1016/j.optcom.2007.11.071
29 Meng X G Wang Z Wang J S Fan H Y 2013 J. Opt. Soc. Am.B 30 1614 https://dx.doi.org/10.1364/JOSAB.30.001614
30 Meng X G Wang J S Gao H C 2016 Int. J. Theor. Phys. 55 3630 https://dx.doi.org/10.1007/s10773-016-2992-0
31 Du C X Meng X G Zhang R Wang J S 2017 Chin. Phys.B 26 120301 https://dx.doi.org/10.1088/1674-1056/26/12/120301
32 Meng X G Wang J S Liang B L Du C X 2018 J. Exp. Theor. Phys. 127 383 https://dx.doi.org/10.1134/S1063776118080113
33 Xu X L Li H Q Fan H Y 2015 Chin. Phys.B 24 070306 https://dx.doi.org/10.1088/1674-1056/24/7/070306
34 Methta C L 1967 Phys. Rev. Lett. 18 752 https://dx.doi.org/10.1103/PhysRevLett.18.752
35 Puri R R 2001 Mathematical Methods of Quantum Optics Berlin Springer-Verlag p. 21
36 Zhang B L Meng X G Wang J S 2012 Chin. Phys.B 21 030304 https://dx.doi.org/10.1088/1674-1056/21/3/030304
37 Meng X G Wang J S Liang B L 2013 Chin. Phys.B 22 030307 https://dx.doi.org/10.1088/1674-1056/22/3/030307
38 Glauber R J 1963 Phys. Rev. 130 2529 https://dx.doi.org/10.1103/PhysRev.130.2529
39 Klauder J R 1963 J. Math. Phys. 4 1005
40 Wang Z Meng X G Fan H Y 2013 J. Phys. A: Math. Theor. 46 135305 https://dx.doi.org/10.1088/1751-8113/46/13/135305
41 Meng X G Xu X J 2014 Int. J. Theor. Phys. 53 1239 https://dx.doi.org/10.1007/s10773-013-1921-8
42 Fan H Y Jiang T F 2007 Mod. Phys. Lett.B 21 475 https://dx.doi.org/10.1142/S0217984907012943
43 Fan H Y Xu X F 2004 Opt. Lett. 29 1048 https://dx.doi.org/10.1364/OL.29.001048
44 Meng X G Goan H S Wang J S Zhang R 2018 Opt. Commun. 411 15 https://dx.doi.org/10.1016/j.optcom.2017.11.005
45 Gardiner C W Zoller P 2000 Quantum Noise 2 Berlin Springer-Verlag p. 147
46 Meng X G Wang J S Liang B L Han C X 2018 Front. Phys. 13 130322 https://dx.doi.org/10.1007/s11467-018-0856-1
47 Xiang S H Zhu X X Song K H 2018 Chin. Phys.B 27 100305 https://dx.doi.org/10.1088/1674-1056/27/10/100305
48 Wang Z Meng X G Li H M Yuan H C 2014 Int. J. Mod. Phys.B 28 1450115 https://dx.doi.org/10.1142/S021797921450115X
49 Wang J S Meng X G Fan H Y 2015 Chin. Phys.B 24 014203 https://dx.doi.org/10.1088/1674-1056/24/1/014203
50 Wang J S Meng X G Fan H Y 2011 Chin. Phys. Lett. 28 104209 https://dx.doi.org/10.1088/0256-307X/28/10/104209
51 Meng Y Wang J S Meng X G 2011 Int. J. Theor. Phys. 50 3348 https://dx.doi.org/10.1007/s10773-011-0838-3
52 Meng X G Liu J M Wang J S Fan H Y 2019 Eur. Phys. J.D 73 32 https://dx.doi.org/10.1140/epjd/e2018-90224-6